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ABSTRACT 

We report the discovery of a microlensing exoplanet OGLE-2012-BLG-0563Lb with the planet-star 
mass ratio of ^ 1 x 10“^. Intensive photometric observations of a high-magnification microlensing 
event allow us to detect a clear signal of the planet. Although no parallax signal is detected in the 
light curve, we instead succeed at detecting the flux from the host star in high-resolution JHK'- 
band images obtained by the Subaru/KOISS and Infrared Camera and Spectrograph instruments, 
allowing us to constrain the absolute physical parameters of the planetary system. With the help 
of spectroscopic information about the source star obtained during the high-magnification state by 
Bensby et al., we find that the lens system is located at 1.3 ® kpc from us, and consists of an M 

dwarf (0.34 ^ 0(20 ^ 0 ) orbited by a Saturn-mass planet (0.39 ^ 0^23 Mjup) at the projected separation 
of 0.74 io ;42 AU (close model) or 4.3 ^ 2^5 AU (wide model). The probability of contamination in the 
host star’s flux, which would reduce the masses by a factor of up to three, is estimated to be 17 %. 
This possibility can be tested by future high-resolution imaging. We also estimate the (J — Ks) and 
[H — Ks) colors of the host star, which are marginally consistent with a low metallicity mid-to-early 
M dwarf, although further observations are required for the metallicity to be conclusive. This is the 
fifth sub-Jupiter-mass (0.2 < rUp/Mjup < 1) microlensing planet around an M dwarf with the mass 
well constrained. The relatively rich harvest of sub-Jupiters around M dwarfs is contrasted with a 
possible paucity of ^1-2 Jupiter-mass planets around the same type of star, which can be explained 
by the planetary formation process in the core-accretion scheme. 

Keywords: planetary systems — planets and satellites: detection — planets and satellites: gaseous 
planets — stars: late-type — techniques: high angular resolution — techniques: photo¬ 
metric 


1. INTRODUGTION 

Microlensing is a unique and powerful technique to 
probe exoplanets with a wide range of masses just be¬ 
yond the snow line, where gas-giant planets can effi- 
ciently form according to the core-accre tion models (e.g., 
iPollack et al.lll99^lKoknbo fc Idall2002tl . The planetary- 
mass distribution probed by microlensing therefore pro¬ 
vides valuable information about the planetary formation 
process, which is less affected by several post-formation 
effects such as orbital migration and mass loss due to 
stellar irradiation. In addition, microlensing is most sen¬ 
sitive to exoplanets around M dwarfs including late-type 
ones, which have not sufficiently been surveyed by other 
detection techniques due to their faintness. The core¬ 
accretion models predict that massive Jovian planets are 
rare around low-mass stars due to the lack of planet¬ 
forming materials fe.g.. llda fc LiDl[2005[l . which can thus 
be tested by microlensing. 

Thanks to a huge effort by microlensing surveys and 
follow-up projects to date, the number of microlensing 
planets has reached 350 among which ^ 60% are hosted 
by M dwarfs. These discoveries have revealed that low- 
mass planets are much more abundant than massive ones, 
whic h is in agreement with the core-accretion scenar¬ 
ios (Gould et al inoia; iSumi et ^ 120101 : iGassan et al.l 
I 2 OI 2 I. On the other hand, super-Jupiter-mass planets 
(> 2 M.T„n) have also been discovered around M dwarfs 

(e.g.. l^ng et al.ll2009l : lBatista et aT]l20lHlTsaDras et al.l 

^ Microlensing Observations in Astrophysics (MOA) Collab¬ 
oration. 

^ Microlensing Follow-up Network (fiFUN) Collaboration. 

^ Optical Gravitational Lensing Experiment (OGLE) Collab¬ 
oration. 

^ Probing Lensing Anomalies NETwork (PLANET) Collabo¬ 
ration. 

^ RoboNet Collaboration. 

® http://exoplanet.eu 


l20I4f) . which challenges the same scenarios. 

However, the statistics of microlensing planets are 
not yet high enough to draw a clear structure of the 
planetary-mass distribution, in terms of number and ac¬ 
curacy. In particular, about half of all planetary mi¬ 
crolensing events do not show parallax effects in the light 
curves, without which one cannot measure the absolute 
masses of the planet and host star from the light curve 
alone. In such cases, the physical parameters of the plan¬ 
etary system have often been estimated by the Bayesian 
technique, which uses Galactic-model priors (a stellar- 
mass function, stellar number density, and stellar veloc¬ 
ity distribution) to draw posterior probability distribu¬ 
tions of the physical parameters. This technique could be 
meaningful if the number of planets is statistically large 
enough, however, the individual values are not accurate. 
Furthermore, this technique relies on an assumption that 
the planet occurrence probability is uniform for all stars 
independent of stellar properties, such as stellar mass 
and Galactic location, and therefore the results should 
be treated with caution. 

Another method to constrain the physical parameters 
of the lens system is detecting (or putting an upper 
limit on) the emission from the host (lens) star by high- 
resolution imaging. High resolution is essential to de¬ 
blend unrelated stars and extract the lens-|-source com¬ 
posite flux. Although, for the current facilities, it is usu¬ 
ally not possible to spatially resolve the lens star from 
the background source star until a decade after the mi¬ 
crolensing event, even without resolving the two stars, 
the lens star’s flux can be extracted by subtracting the 
source star’s flux (obtained by a light-curve analysis) 
from the lens-l-source composite flux. The extracted lens 
flux provides a mass-distance relation of the host star, 
allowing us to solve for the mass and distance by com¬ 
bining with another mass-distance relation provided by 
the angular Einstein radius 9e, which can be derived in 
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most planetary microlensing events. 

High-resolution imaging has now become an important 
tool to constrain the physical parameters of the lens sys¬ 
tems. So far, 11 planetary events have been imaged at 
high resolution usin g the Hubble Space Telescope (e.g., 
iBennett et al.l 1200611 or ground-based adaptive-optics 
(AO) instruments ( Very Large Telescope: VLT/NACO 
or Keck/N1RC2; e.g.. lJanczak et al]l2010l : [Bennett et al.l 
[Moll . Among them, five events did not show par¬ 
allax effects in their light curves, and therefore the 
high-resolution images played a complementary role to 
constrain the lens-system parameters. Regarding M- 
dwarf host star events, the following six have been im¬ 
aged to_dnte_wd^ high resolutiom^OGLE-2005-BLG- 
071 (lUdalski et aLllMalPoM et al.ll2009ll OGL E-2006- 
BLG-109 (iGaudi e t al.ll200^ Bennett e t al.ll2010ll MOA - 
2007-BLG-192 ([Bennett et al.l 200^ Kubas et al.l 1201211 , 
and MOA-2009- BLG-387 (iBatista et al.l 1201 111 . OGLE- 
2013-BLG-0341 (iGould et al.ll20l4 Bati sta et al. in 
prepa ration), and MOA-2011-BLG-262 (IBennett et al.l 
1201411 . Among them, the first five events also showed 
parallax signals in the light curves, and thus the high- 
resolution images have been used to support or reinforce 
the results from the light-curve analyses. As for the last 
event, MOA-2011-BLG-262, the lens flux was not clearly 
detected. This event has actually two degenerate solu¬ 
tions: a planet-M-dwarf system in the Galactic bulge and 
a nearby free-floating planet orbited by a moon. 

In this paper, we report the discovery of a new plane¬ 
tary microlensing event with the planet-star mass ratio 
of ~ 10“^ without parallax signal, for which we constrain 
the physical parameters of the lens system by combining 
light-curve analysis with Subaru/AO imaging. This is the 
first M-dwarf-host planetary event without parallax for 
which the lens flux is clearly detected, and therefore the 
AO imaging plays a crucial role for deriving the physical 
parameters. 

The rest of this paper is organized as follows. The ob¬ 
servations and light-curve analysis are described in Sec¬ 
tions [2] and El respectively. The properties of the source 
star are investigated in Section jT] Sections |S] and [S] de¬ 
scribe the extraction of excess flux on the source star 
and constraint on the physical parameters of the lens 
system, respectively. We discuss the results in Section [7] 
and summarize in Section 8. 


2. OBSERVATIONS 
2.1. The Microlensing Event 

The microlensing event OGLE-2012-BLG-0563 oc¬ 
curred on a star located at the equatorial coordinate of 
(a, ^)j 2 ooo = (18^05“57".72, -27“42'43"2) and the Galac¬ 
tic coordinate of (/, b) = (3?31, -3?25). This event was 
first discovered by the OGLE collaboration on 2012 May 
1 UT (JD' = JD — 2450000 = 6049), during its reg¬ 
ular photometric monitoring toward the Galactic bulge 
by using the 1.3m Warsaw telescope equipped with the 
wide field (1.4deg^) ca mera OGLE-IV at L as Campanas 
Observatory in Ghile ([Udalski et al.l 1^151) . The event 
was discovered in the field BLG519, which was regu¬ 
larly monitored twice a night with an /-band filter. The 
same event was independently discovered as MOA-2012- 
BLG-288 by the MOA collaboration on 2012 May 18 UT 
(JD^ = 6066), by using the 1.8m MOA-II telescope and 


the w ide field (2.2-deg^) camera moa-cam3 (|Sako et al.l 
120081) at Mt. John University Observatory (MJUO) in 
New Zealand. MOA monitored the event (in field gbl4) 
with a typical cadence of 20 minutes using a custom R+I 
filter. After these discoveries, both teams have made the 
light curves public and updated them frequently. 

On May 20 (JD' = 6067.9), the /iFUN collaboration 
circulated an alert that the event was peaking at a high 
magnification, which means that the source star was ap¬ 
proaching very close to the lens star on the sky plane. 
In such a case, there is a high probability that the light 
curve will show anomalous features around the peak if 
the lens star hosts a planet, because one of the caus¬ 
tics (central caustic) produced by a planetary system 
is always created near the host star on the sky plane 
(iGriest fc Safizad^ 119981 : iRhie et all 120001 : |H an fc Kiinl 

[2001)7 ^ 

After the circulation of the high-magnification alert, 
the /xFUN and RoboNet collaborations started to follow 
up the event with high cadence photometry. //FUN used 
the following telescopes (filters): the 1.3m GTIO tele¬ 
scope in Ghile {V, I, and H), the 0.40m Auckland tele¬ 
scope in New Zealand (/?), the 0.36m Possum telescope 
in New Zealand (/?), the 0.36m Farm Gove telescope in 
New Zealand (clear filter), and the 0.30m PEST in Aus¬ 
tralia (clear filter). The 1.3m GTIO’s //-band data were 
simultaneously obtained with the V- or /-band data with 
the same instrument. RoboNet used three 2.0m tele¬ 
scopes, including the Liverpool Telescope (LT) in Ca¬ 
nary Islands, Spain (/'); the Faulkes Telescope North 
(FTN) in Hawaii, USA (*'); and the Faulkes Telescope 
South (FTS) in Australia (/'). In addition, OGLE and 
MOA increased the observational cadence with their sur¬ 
vey telescopes. MOA also utilized the 0.61m B&C tele¬ 
scope at MJUO with /-band hlter to follow up the event. 
On May 21 13:30 UT (JD' = 6069.06), MOA reported 
the detection of an anomaly as a result of its high ca¬ 
dence observations. The event peaked at around May 21 
13:20 (JD' = 6069.06) with the maximum magnification 
of more than 600, around which the event was well cov¬ 
ered by several telescopes in New Zealand and Australia. 
As a result, a quick look light curve clearly showed an 
anomalous asymmetric feature around the peak, meaning 
that the lens star hosted a companion. Prompt analyses 
of the light curve by several teams indicated that this 
was a firm planetary event having the planet-star mass 
ratio of ~ 10“^. We summarize these photometric obser¬ 
vations in Table [1] and show the observed light curve in 
Figured! We note that data obtained by Possum, Farm 
Cove, and FTN are not used for the analysis due to the 
following reasons; data from Farm Cove do not constrain 
the model because they consist of just two short epochs, 
and data from Possum and FTN are of relatively low 
quality. 

In addition to the light curve, a high-re solution spec¬ 
trum of the source star was obtained by iBensbv et al.l 
([Mil with VLT/UVES on 2012 May 19 (JD' = 6067), 
when the source was magnified by a factor of ~50. They 
reported that the source star was a metal-poor G dwarf 
with Tes = 5907 ± 89 K, logg = 4.40 ± 0.10, and 
[Fe/H]= -0.66 ±0.07. 

After the event was over, all the observed photometric 
data were carefully re-reduced by their respective teams. 
For the re-reductions, OGLE and MOA used their cus- 
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Figure 1. (Top) light curve of the event OGLE-2012-BLG-0563. 
The correlation between color and instrument is shown on the right. 
The solid black line indicates the best-fit planetary microlensing 
model. The inset shows the geometry of the event projected onto 
the lens plane. The origin is the centroid of the planetary sys¬ 
tem, and the X axis is defined as the host star-planet axis. The 
black and gray lines indicate the central caustic for the best-fit 
close model (which is almost identical to that for the wide model) 
and the source trajectory, respectively. The blue circle represents 
the source star at to, with the size showing the source size. The 
black dots are tick marks for time, with number indicating HJD- 
2456068. (Middle) a zoom around the peak. (Bottom) residuals of 
the zoomed light curve from the best-fit model. 


tomi zed pipeline s as d escribed in lUdalski et al.l (|2015[ ) 
and iBond et ^ (120011). respectively; /x FUN used the 
DoPHOT package (|Schechter et al.lll993D for the CTIO 
1.3m iJ-band data and a variant of the PySlS package 
(|Albrow et al.l[200^ for the others; a nd RoboNET used 
the DanDIA package (jBramichl 120081 ). All the pipelines 
except for DoPHOT apply the Difference Image Anal¬ 
ysis (DIA) method for photometry in order to achieve 
good photometric precisions in the star-crowding field. 
Note that because a bright neighboring star was located 
2" 2 away from the target star, the photometry was care¬ 
fully done to minimize systematics (e.g., by masking the 
bright star). The time for each data point was converted 
to heliocentric Julian day (HJD) based on UT. 


star having i?=ll and being separated by 26" from the 
source star. The field was observed through J-, 77-, and 
T^'-band filters, each with 30 s exposure x 15 times at 
five dithering points within 2" (100 pixel) square. The 
airmass toward the target field was 1.74-1.54 during the 
observations. The natural seeing was ^0"5-0"6, and the 
AO-worked seeing was 0"23-0'.'30, 0"19-0'.'24, and 0"17- 
0725 for J, 77, and K', respectively. 

All the J-, 77-, and 7T'-band AO images were dark- 
subtracted and flat-fielded in a standard manner. For 
the flat fielding, we used sky flat images obtained during 
evening twilight on the observation night. All the images 
in each band were combined in average, after the image 
positions were aligned and sky levels were subtracted. 
The image-overlapping region was trimmed for further 
analyses, which resulted in the effective FOV of 18" x 
18". 


3. LIGHT CURVE ANALYSIS 

3.1. Model Description 

Because the light curve of the event clearly shows an 
asymmetric feature around the peak that cannot be ex¬ 
plained by a standard single-lens microlensing model, 
we introduce a binary-lens microlensing model assum¬ 
ing that the lens consists of two objects. The standard 
binary-lens model can be described with seven basic pa¬ 
rameters: the time of the closest source approach to the 
binary centroid, fo; the Einstein radius (7?e) crossing 
time, t-E] the minimum impact parameter, uq] the mass 
ratio of the binary components, q; the projected sepa¬ 
ration of the binary components in units of 7?e, s; the 
angle between the source star trajectory and the binary- 
lens axis, a; and the angular source radius in units of 
6 *e, P- In addition, the model requires two instrument- 
dependent parameters of the unmagnified source flux 
and the blended flux 7 b- We model the intensity dis¬ 
tribution on the surface of the source star with a linear 
limb-darkening raw of 7(0) = 7(0) [1 — rtjc(l — cos0)], 
where 9 is the angle between the normal to the stellar 
surface and the line of sight, 7(0) is the stellar intensity 
as a function of 0, and ux is a coefficient for filter X. 
For My . we adopt the theoretical values of iClaret et al.l 
(j2013f ) for a G dwarf with the temperature of 5900 K and 
log (?=4.5, namely, uv=0.656, uh=0.572, u^+/= 0.528, 
Uj,=0.504, u/=0.483, and u//=0.290. The ur+i value is 
calculated as the mean of ur and uj. For clear filter, 
we adopt ur as an approximation. The model calcu¬ 
lation is done by using a customized code developed by 
MOA, w hich is based on the image-centered ra y-shooting 
method (|Bennett fc Rhidri996l : lBennetfll2010D . 


2.2. Subaru AO Imaging 

We conducted high-resolution imaging of the event 
field by using the 8.2m Subaru telescope equipped with 
the AO instrument A0188 and Infrared Gamera and 
Spectrograph flRCS. lKobavashi et al.l[2000l) at 6:06-7:03 
UT on 2012 July 28 (JD' = 6137.8), when the source star 
was still magnified by a factor of 1.47. We used the “high- 
resolution” mode of IRCS, which provides a pixel scale 
of 20.6mas pixeU^ and a field of view (FOV) of 2171 
X 2171. We used the natural guide star (NGS) mode 
for AO. As an NGS, we selected a bright neighboring 


3.2. Error Normalization 

The initially calculated flux uncertainties in the light 
curves are normalized as follows. First, all the light 
curves are simultaneously fit with the binary-lens mod¬ 
els by th e Markov Ghain M onte Garlo (MCMC) method 
following iSumi et al.l (j20ir)D . Then the flux uncertainties 
in each light curve are scaled by a single factor so that 
the per degrees of freedom, for each light curve 
becomes unity. This process is iterated until the best-fit 
light curve model becomes converged. Next, after exclud¬ 
ing 4(7 outliers, we renormalize the flux uncertainties for 
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Table 1 

Summary of Observations 


Telescope 

Diameter (m) 

Filter 

A^use ’‘/Mobs 

The OGLE collaboration 

OGLE 

1.3 

I 

700/751 

The MOA collaboration 

MOA-II 

1.8 

R + I 

5449/5947 

B&C 

0.61 

I 

385/439 

The uFUN collaboration 

CTIO 1.3m 

1.3 

V 

12/12 

CTIO 1.3m 

1.3 

I 

30/31 

CTIO 1.3m 

1.3 

H 

110/140 

PEST 

0.30 

clear 

163/172 

Auckland 

0.40 

R 

44/46 

Possum 

0.36 

R 

0/21 

Farm Cove 

0.36 

clear 

0/23 

The RoboNET collaboration 

FTS 

2.0 

i' 

195/218 

LT 

2.0 

i' 

20/23 

FTN 

2.0 

i' 

0/10 


^ The number of data points used in the analysis. 
^ The number of observed data points. 


each light curve using the following formula, 

+ eLn> (1) 

where ai is the initial error bar of the *th data point in 
magnitude, and k and emin are coefficients for each data 
set. Here, the term Cmin represents systematic errors 
that dominate when the stellar brightness significantly 
increases. We adjust the k and emin values so that the 
cumulative distribution sorted by magnitude is close 
to linear, and Xred becomes unity. Finally, all the normal¬ 
ized light curves are fit again and the flux uncertainties 
are rescaled by a single factor so that data. 

becomes unity. 

3.3. Best-fit Models 

In modeling binary microlensing light curves, several 
models are often degenerate. In particular for central- 
caustic crossing and approaching events, a severe de¬ 
generacy between s and often occurs because the 
central-caustic pattern created by a binary lens with s 
and t hat with ca n be very similar (wide-close degen¬ 
eracy; |Domini3[T99^. To check for this type of degener¬ 
acy and possible other degeneracies, we calculate an 
map in the logq and logs plane by dividing logg ([-5, 0]) 
and logs ([-0.7, 0.7]) into 50 x 40 grids and fitting the 
light curve while fixing log q and log s at each section. 
These logg and logs ranges are chosen such that the 
usual sensitivity region of microlensing is well covered. 
The calculated x^ map is shown in Figure [21 The map 
shows that there are two local minima at (logs, logg) 
^ (-0.4, -3) and (0.4, -3), indicating that the wide-close 
degeneracy clearly exists. On the other hand, there is 
no other local minimum in the map, and the g value at 
the two local minima of ~ 10“^ is well within the plane¬ 
tary range (g < 10“^ for G- and later-type dwarf hosts) 
with Ax^ > 1000 over the models of g > 10“^. There¬ 
fore, the planetary nature of this event is quite robust. 
Note that one can wonder whether there are any other 
degenerated binary (large-mass ratio) models outside the 
searched log s range. However, a very close binary with 


log s < —0.7 would cause a higher-order effect in the light 
curve, due to its orbital motion, which we do not detect. 
In addition, we also check for very wide binary models 
with 0.7 < logs < 2, but find no comparable models 
with the planetary one. 

To properly derive the best-fit model parameters and 
their uncertainties for the wide and close models, we re¬ 
run the MCMC program by starting with the parameter 
values at each local minimum and letting all parameters 
be free. For each model, we conduct 40 independent 
MCMC runs with 3200 steps each, and create poste¬ 
rior probability distributions of the parameters from the 
merged ^10^ steps. The best-fit value and its Icr lower 
(upper) uncertainties for each parameter are calculated 
as the median and 15.9 percentile (84.1 percentile) val¬ 
ues of the posterior probability distribution, respectively. 
The best-fit light curve and caustic models are shown in 
Figure (U while the resultant parameter values, their un¬ 
certainties, and the minimum-x^ values for the wide and 
close models are listed in Table [2] The difference of the 
minimum-x^ values of the two models is only 0.6, mean¬ 
ing that these two models are indistinguishable. Note 
that all the parameter values except for s are almost 
identical between the two models, and we proceed with 
further analyses of only the close model as representative 
of the two, unless s is relevant. 

3.4. Searching for the Microlens Parallax Effect 

In order to know the absolute masses and projected 
separation of the planetary system instead of g and s, 
one needs to derive the total mass of the planetary sys¬ 
tem, mL, and the distance from the Earth to the system, 
Di^. One way to do this is to detect the microlens par¬ 
allax effect in the light curve, which is seen as a slightly 
asymmetric distortion over the non-parallax light curve. 
The distortion depends on ttil, the vectors of the 
Earth’s orbital acceleration, and the lens-source relative 
proper motion ([Gouldl [19921) . From this effect, one can 
obtain an additional parameter tte, which is defined as 
the ratio of 1 AU to the proje cted Einstein radius onto 
the Earth’s plane (|Gouldll2000r) . From tte and the angu- 
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Table 2 

Results of the MCMC Analysis. 


Parameter 

Units 

Close 

Wide 

to 

HJD-2450000 

6069.02790 ± 0.00026 

6069.02811 ± 0.00027 


days 

77.5 ± 2.2 

77.7 ± 2.1 

Uq 

10-3 

1.405 ± 0.040 

1.403 ± 0.039 

q 

10-3 

1.086 ± 0.040 

1.085 ± 0.039 

S 

0E 

0.4134 ± 0.0032 

2.425 ± 0.019 

a 

rad 

5.7843 ± 0.0013 

5.7841 ± 0.0013 

p 

lO-'^ Be 

o 97 +0.24 
-0.26 

q 97 +0.24 
-0.26 

■ 

^min 


7068.6 

7069.2 



Figure 2. A map in the logs vs. logq- plane. 


lar Einstein radius, which will be measured in Section 
m one can derive tol and Dl as 


mL = ■ 


ktte 


nL = 


AU 


TTEfi^E + 1/ns ’ 


( 2 ) 

(3) 


where k = AGjc^ ~ 8.144 mas Mg^ and Ds is the source 
distance. Here G is the gravitational constant and c is 
the speed of light. 

To search for the microlens parallax effect, we fit the 
light curve by freeing two additional parameters of tte^e 
and tte.n, which are the east and north components of 
tte, respectively, where tte is a vector whose length is tte 
and direction is the same as the lens direction relative to 
the source. 

As a result, the parallax model fit gives the best-fit val¬ 
ues of tte.e =-0.024 ± 0.058 and 7rE.N=0.62 ± 0.17 with 
the of 7054.1, which has the improvement of 14.5 com¬ 
pared with that for the non-parallax model. Statistically 
speaking, this is a marginal (S.Scr) detection of the par¬ 
allax signal, however, we consider it suspect because of 
the following reasons. 

The difference between the best-fit parallax and non¬ 
parallax models emerges around the wings of the light 
curve, where the MOA-II and OGLE data are dominant. 
In Figure[31 we show the MOA-II and OGLE light curves 
along with the best-fit parallax (green) and non-parallax 
(cyan) models in the top panel, the residuals of the ob¬ 


served light curves against the non-parallax model in the 
middle panel, and the difference of cumulative be¬ 
tween the two models in the bottom panel. The resid¬ 
ual plot shows that the difference of the two models is 
quite small compared to the error bars of the data points. 
In such a case, systematics in a small number of data 
points can lead to a false positive detection. In fact, 
as indicated by the Acumulative-x^ plot, the major 
improvement comes from the MOA-II data in the lim¬ 
ited range of 6070 < HJD-2450000 < 6110, although the 
model difference emerges more widely over several hun¬ 
dred days. Moreover, the Acumulative-y^ plots for the 
MOA-II and OGLE data are anti-correlated rather than 
correlated, indicating that the detected parallax signal is 
suspicious. Therefore, the detected signal is probably a 
false positive, and we only calculate an upper limit on 
tte < 1.53 (3 a). We note that other effects that can 
mimic the parallax signal, such as the orbital motion of 
the source and/or lens systems, can also be rejected for 
the same reasons. 

4. PROPERTIES OF THE SOURCE STAR 

The 0E value can be derived by dividing p by the angu¬ 
lar radius of the source star 0*, which can be estimated 
from the intrinsic color and magnitude of the source star. 
For most microlensing events, no spectroscopic informa¬ 
tion of the source star is available, due to the intrinsic 
faintness of the star. In this case, the intrinsic color and 
magnitude of the source star are estimated by using red 
clump stars in the Galactic bulge as a standard candle, 
assuming that the dust extinction for the source star is 
the same as that for those red clump stars. We first 
derive 0* using this “standard” method in Section 14.11 
On the other hand, in the current case, a high-resolution 
spectrum of the source star was obt ained during the high - 
magnification state of the event by iBensbv et al.l (1201311 , 
as part of a systematic elemental-abundance study of the 
bulge stars. With this spectral information, the source’s 
intrinsic color and magnitude can be derived with a mini¬ 
mum assumption about the extinction for the source star. 
We show how the spectral information improves the 9^ 
estimation in Section 14.21 In addition, we measure the 
distance to the source star by combining the photometric 
and spectroscopic information in Section 14.31 

4.1. Angular-radius Estimation Using a Standard 
Method 

We first derive the apparent I- and U-band magni¬ 
tudes of the source star, Ig and Vs: respectively, from 
the light-curve fitting. We derive Is from the two light- 
curve data sets obtained by the OGLE and CTIO 1.3m 
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Figure 3. (Top) The OGLE (black) and MOA-II (red) light 
curves along with the best-fit models with (green) and without 
(cyan) the parallax effect. (Middle) Residuals from the non¬ 
parallax best-fit model. The data points are binned in 10 days 
for clarity. (Bottom) Difference of the cumulative between the 
parallax and non-parallax models for the OGLE (black) and MOA- 
II (red) data. A negative value means a preference for the parallax 
model. Note that, although this plot indicates that the total 
improvement from these two data sets is ^^38, that from all data 
sets is 14.5, as described in text. 


telescopes. Calibrating the instrumental magnitude to 
the standard (Landolt ) one via the OGLE-III catalog 
(|Szvmahski et al.] 1201 111 , we obtain Is = 20.147 ± 0.031 
and Is = 20.082±0.031 from the respective data sets. By 
taking the mean of the two values, we obtain a final value 
of Is = 20.115 0.031. Note that the uncertainty of Is 
is conserved in this calculation because this uncertainty 
is dominated by the light-curve model uncertainty and 
therefore the two Is values are correlated. Vs is derived 
from the V-band light curve obtained by the CTIO 1.3m 
telescope. In the same way as for the 1 band, we derive 
the calibrated source magnitude of Vs = 21.595 ± 0.032. 
Then, from Is and Vs, we obtain the apparent source 
color of {V - I)s = 1.480 ± 0.032. 

Next, we estimate the extinction and reddening for the 
source star using red clump stars in the Galactic bulge. 
In FigurelH we plot the color magnitude diagram (CMD) 
toward the event coordinate ('^218 x 218 area) created 
from the OGLE-111 catalog, along with the measured Is 
and (V — I)s- The source position is largely consistent 
with, but a bit fainter than, the region of MS stars in the 
Galactic bulge, which is consistent with the fact that the 
source is located behind the cen ter of the bulge along 
the line of sight (see Section 1431) . The centroid of the 
red clump stars on the CMD is measured as (V — /)rc = 
1.883 ± 0.013 and Irc = 15.361 ± 0.015, respectively. 
These values are then compared to the intrinsic color 
and magnitu de of the red clump stars of (E — /)rc o = 
1.06 ± 0.06 (IBensbv et all [Mlh and /rc,o = 14.34 ± 


14 

1 6 

18 

20 

1 1.5 2 2.5 

(V-I) 


-|-Red clump 

• Source 



Figure 4. Color-magnitude diagram (CMD) toward the event 
field (~ 2^8 X 2.^8 area) created from the OGLE-III catalog. The 
centroid of red clump giants and the source position are indicated 
by the red cross and blue circle, respectively. 


0.04 (jNataf et al.l l2013ll . leading to the reddening and 
extinction magnitudes for the red clump stars oi E{V — 
I) = 0.82 ± 0.06 and Aj = 1.02 ± 0.04, respectively. 
Assuming that this reddening and extinction can also 
be applied for the source star, we derive the intrinsic 
color and magnitude of the source star as (V — /)s,o = 
0.66 ± 0.07 and Isfi = 19.10 ± 0.05, respectively. 

Finally, we estimate the angular radius of the source 
star 0* from /s,o and (E —/)s,o using the following equa¬ 
tion: 


log Old = 0.5014 -f 0.4197(E -I)- 0.2/, (4) 

where 0 ld is the limb-darkened stellar angular diame¬ 
ter. This linear equation is derived from a subset of 
th e interferometri c ally m easured stellar radii presented 
in iBovaiian et all (|2014f) . restricting stars with 3900 K 
< Teff < 7000 K to improve the ht for FGK stars. This 
gives 0* = 0 ld/ 2 = 0.454 ± 0.033 ^as. 

The derived source properties are summarized in Table 

El 


4.2. Angular-radius Estimation with Spectral 
Information 

4.2.1. From {V — I) and I 

The spectroscopically measured effective tempera¬ 
ture Teff = 5 907 ± 89 K and rn etallicity [Fe/H] = - 
0.66 ± 0.07 (|Bensbv et al.l 1201311 of the source star 
can be directly converted to the source’s intrinsic 
col or. Using the col o r-met allicity-temperature relation 
of ICasagrande et al.1 (|2010ll . we derive (E — /)s,o = 
0.660 ± 0.027. This value is in good agreement with the 
photometrically derived (E — /)s,o but has a smaller un¬ 
certainty by a factor of ^2.5. Subtracting this value from 
the apparent source color of (E — I)s = 1.480 ± 0.032, 
we derive the reddening magnitude up to the source star 
of E{V — I) = 0.820 ± 0.042. Then, assuming that the 




































Table 3 

Properties of the Source Star 


Parameter 


Value 

Value 



(w/o spec, info) 

(w/ spec, info) 

Teff (K) 
[Fe/H] 
logs (cgs) 


— 

5907 ± 89 
-0.66 ± 0.07 

4.40 ± 0.10 

V's 


21.595 ± 0.032 

— 

-^s 


20.115 ± 0.031 

— 

Hs 


18.569 ± 0.032 

— 

(V" - Ds,o 
iV-H)s,o 
Fs,o 


0.66 ± 0.07 

19.76 ± 0.08 

0.660 ± 0.027 

1.41 ± 0.06 
19.73 ± 0.07 

fs,o 


19.10 ± 0.05 

19.14 ± 0.08 

Hs.o 


— 

18.314 ± 0.041 

F(V-l) 


0.82 ± 0.06 

0.820 ± 0.042 

E(J-Kg) 

Av 


0.244 ± 0.109 ^ 
1.84 ± 0.07 

0.266 ± 0.019 
1.80 ± 0.09 

Ai 


1.02 ± 0.04 

0.98 ± 0.08 

Aj 


— 

0.423 ± 0.024 

Ah 


— 

0.255 ± 0.029 

Ak, 


— 

0.157 ± 0.017 

Ot [pas] [1, V — 
e, [pas] {H, V - 

U 

H) 

0.454 ± 0.033 

0.446 ± 0.023 

0.444 ± 0.014 

Mv 

Mi 

Mh 
t (Gyr) 

Ds (kpc) 


— 

-0.14 

4 28 

_o.i9 

3.54 

13 

Q 1 +0-9 
-1.1 

® From IBensbv 

et 

all (1201311. 



The value adopted for the rest of analyses. 


4.2.2. From {V - H) and H 

We recalculate 0* from {V — H) and H of the source 
star in order to check for the robustness of the previous 
estimation. In general, one can obtain a better estima¬ 
tion of 0* from (V — H) and H rather than from {V — I) 
and I (c.f., Table 4 of iKervella et al.|l2004f) . 

In the same way as for {V — d)s, 0 i we calculate the in¬ 
trinsic {V — H) source color of {V — i?)s,o = 1-41 ± 0.06, 
where H is in the 2MASS system. Then, the intrin¬ 
sic i?-band source magnitude i?s,o is estimated as fol¬ 
lows. First, the apparent i/-band source magnitude is 
measured from the microlensing fit to the i7-band light 
curve obtained by the CTIO 1.3m telescope, yielding 
Hs = 18.569 ± 0.032 in the 2MASS system. Next, the 
Fl-band extinction Ah is estimated by combining the 
(J — Ks) reddening of E{J — Kg) = Rjkvi x A(14 — 
I) = 0.266 ± 0.019 (adopting O.Ol for the uncertainty 
of Rjkvi) and two extinction coefficients of AjjAx^ 
and Ah/Ak^- We adopt Aj/Ak^ = 2.70 ± 0.15 and 
Ah/Ak^ = 1.63 ± 0.05, which ye the mea n valu es for 
the Galactic bulge estimated bv iChen et al.l (|2013ll . but 
have conservative uncertainties taking into account the 
non-uniformity of these coefficients toward the Galac¬ 
tic bulge. We solve these equations for Ah and de¬ 
rive Ah = 0.255 ± 0.029. As a by-product, we also 
obtain Aj = 0.423 ± 0.024 and Ak, = 0.157 ± 0.017, 
which will be used in Section 15.31 Finally, we derive 
Hs.o = Hs -Ah = 18.314 ± 0.041. 

The derived {V — H)s.q and i7s,o values are then con¬ 
verted to 0* using the following equation: 


following extinction law 
Ai = 1.217 X 

E{V - /)[1 -k 1.126 X {Rjkvi - 0.3433)] (5) 

(|Nataf et al.l l2013ll applies along the line of sight, and 
adopting Rjkvi = E{J - Kg)/E(V - I) = 0.3249 for 
the e vent coordinate fr om the online Extinction Calcula- 
toiQ (|Nataf et al.ll2013[l . we obtain Aj = 0.98 ± 0.08 (the 
error includes the uncertainty of Equation ([5]) , for which 
we adopt 0.06 mag). Consequently, we derive the intrin¬ 
sic source magnitude of /s,o = Is — Aj = 19.14 ± 0.08. 
This value is also consistent with the previous estimation 
but has a bit larger uncertainty, which mainly comes from 
the uncertainty of the estimation of Aj. Note, however, 
that this Aj estimation is not based on any assumption 
about the absolute extinction or reddening, but on the 
extinction law of Equation ([5]). The source angular ra¬ 
dius is then calculated by using Equation (|4]) as 

0* = 0.446 ± 0.023 /ias, (6) 

whose uncertainty is 1.4 times smaller than the previous 
one. Using this 0* and p, we derive the angular Einstein 
radius of 

Be = l-36lo;]2 mas. (7) 

We also derive the geocentric source-lens relative proper 
motion from Be and Ie as 

Pgeo = 6.4to;5 mas yr“\ (8) 


^ http://ogle.astrouw.edu.pl 


log0LD = 0.53598-h0.07427(U- i7) 

-h0.04511[Fe/H]-0.2i7, (9) 

where H is in the Johnson magnitude system. This equa¬ 
tion is derived in the same way as for Equation ([4]), but 
includes the metallicity term because a small metallicity 
dependence appear s in this relation. We adopt [Fe/H]=- 
0.66 ± 0.07 from iBensbv et al.l (|2013ll . We convert 
iy — H) in the 2MASS magnitude to that in the John¬ 
son system via the transformation of {V — Tfjohnson) = 
0.28302xl0-2-hl.0021(U-772MASs)+0.36618xl0-2(V- 
if2MASs)" - 0.17906 X 10-2(U - i72MASs)" + 0.16113 X 
10“^(U—i72MASs)"^i whic h is derived by com bining Equa¬ 
tion (B5) and (B8) from iCarpented (l2001h and a color- 
color relation of Table A5 from iKenvon fc HartmanrJ 
(|1995f) . As a result, we derive 0* = 0.444 ± 0.014 pas, 
which is quite consistent with the previous estimations, 
and even has the smallest uncertainty. However, this 
value is derived partly by relying on the two extinc¬ 
tion coefficients of Aj/A k^ and A h/Ak.. which cany ary 
depending on the line of sight (|Chen et al.l 1201311 al¬ 
though we adopt the mean values for the galactic bulge. 
That could produce some systematics that are difficult 
to asses. Therefore, we adopt the 0* value derived in the 
previous section for the rest of analyses to be conserva¬ 
tive. 

4.3. Distance to the Source Star 

The distance to the source star Dg can be measured 
by estimating the absolute magnitude of the source star 
in addition to the intrinsic one. To this e nd, we use the 
isochrone models of PARSEC version 1.2S l|Bressan et al.l 
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l2012fl . Through the web interface CMD 2.713 we obtain 
a grid of isochrone tables for the stellar age in the range 
of 8.0 < logt[yr] < 10.13 with a step size of logt=0.02, 
and for the metallicity in the range of -0.88 < [M/H] < 
-0.45 with a step size of 0.035. For each table, we further 
grid the table interpolatively by logTes with a step size 
of 0.001. Then, for each grid, we calculate the following 
value 


2 (TeS - Tf 

^ 0-2 

Teff 


([Fe/H] - M)- 

^[Fe/U] 


^logS 


where T, M, and G are the model temperature, metal¬ 
licity, and surface gravity, respectively. We then find 
the best-ht values of the absolute V, /, and H mag¬ 
nitudes My, M/, and Mh^ as well as the stellar age 
t, by minimizing the value. In addition, Itr uncer¬ 
tainties of these parameters are calculated by searching 
for the region where = 1. The resultant values 

and uncertainties are My = 4.991 q 24 , = 4.281q 4 |, 

Mh = 3.54l n4^, a nd t = 13^5 Gyr. We note that 
IBensbv et al.l (|2013ll also estimated My, Mj and t by us¬ 
ing the Yonsei-Yale isochrone models as 4.89, 4.23, and 
10.2^45 Gyr, respectively. Our estimations are consis¬ 
tent with theirs. 

Gombining these absolute magnitudes with the intrin¬ 
sic ones derived in Section we calculate the source 
distance to be Us = 8.8 9.4 tyi, and 9.0 tyo kpc 

for V, I, and H, respectively. By taking the mean of the 
three values while conservatively keeping the uncertain¬ 
ties, we obtain a final value of Dg =9.1 I'j'’® ^pc. On 
the other hand, the distance to the centroid of the red 
clump stars toward the event co ordinate is e s timat ed us¬ 
ing the Extinction Calculator of iNataf et al.l (l2013f) to be 
Urc = 7.9^Q g kpc, which indicates that the source star 
is likely located at the far side of the Galactic bulge. This 
result is in agreement with the fact that the microlensing 
event rate is higher for far-side source stars because of 
the higher stellar density per unit solid angle. 


5. EXTRACTING EXCESS ELUX 

Although the orbital parallax is not detected in the mi¬ 
crolensing light curves, the absolute lens mass and dis¬ 
tance can be constrained if the flux from the lens star 
is detected. Although the lens and source stars cannot 
spatially be resolved at present, the lens flux can be ob¬ 
served as an excess of flux at the source position. To 
do this, we analyze the high-resolution images obtained 
with Subaru/IRCS. 


5.1. Identifying the microlensing target 

The first step to extract the excess flux is to identify 
the microlensing target (source-l-lens superposition) on 
the 5M6arM/IRCS images. To do so, we use an OGLE- 
IV /-band image obtained near the brightness peak with 
the source magnification of more than 400 times (upper 
right panel in Figure [5]), which overwhelms any blended 
fluxes so that the source position can unambiguously be 
measured. Because the target position is almost identical 
to the source position at the time of the Subaru obser¬ 
vation, the measured source position on the OGLE peak 


® http://stev.oapd.inaf.it/cgi-bin/cmd_2.7 



Figure 5. (Upper left) the OGLE 7-band reference image ob¬ 
tained well before the microlensing event. The location of the 
source star to be magnified is indicated by an arrow. (Upper right) 
an OGLE 7-band image obtained near the event peak. (Lower left) 
the IRCS J-band image. (Lower right) a zoom of the image in the 
lower-left panel. The calibrated source position is indicated by cyan 
lines. The region in which the sky value is calculated is shown in 
yellow box. 


image can be used to identify the target on the IRCS 
image once the instrumental coordinate of the OGLE 
image is calibrated to that of the IRCS one. The coor¬ 
dinate calibration is done by constructing a calibration 
ladder via an OGLE reference image (upper left panel 
in Eigure [5]) obtained under a better seeing condition 
(I'.'O) compared to the OGLE peak image (I"2). Specif¬ 
ically, we first calibrate the instrumental coordinate of 
the OGLE peak image to that of the OGLE reference 
one using the centroid positions of 18 bright common 
stars (first calibration), and then calibrate the instru¬ 
mental coordinate of the OGLE reference image to that 
of the IRCS/J-band one using the centroid positions of 
nine well-isolated common stars (second calibration). All 
the stellar centroids are measured by using the DoPHOT 
package. The coordinate calibrations are done by using 
the IRAFI3 geo map algorithm with free parameters of 
xy shifts and rotation for the first calibration, and xy 
shifts, xy pixel-scale magnifications, and rotation for the 
second calibration. The rms values for the first and sec¬ 
ond calibrations are 38 mas for 18 stars and 44 mas for 
nine stars, respectively, meaning that the total Icr cali¬ 
bration error is 17 mas. We measure the source centroid 
position on the OGLE peak image with an adopted Icr 
uncertainty of 0.1 pixel, or 26 mas, which is then cali¬ 
brated to the coordinate on the IRCS image (marked as 
cyan lines in Eigure [5]) with the total uncertainty of 31 
mas. On the IRCS image, we find that there is one stellar 
object close to the source position with a separation of 

® IRAF is distributed by the National Optical Astronomy Ob¬ 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the US National Science Foundation. 
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only 44 mas, which is comparable with the measurement 
uncertainty. We therefore conclude that this object is 
the microlensing target. 

5.2. Photometry of the microlensing target 

The next step is to measure the brightness of the mi¬ 
crolensing target. We perform aperture photometry of 
the target object on the IRCS J-, H-, and T^'-band im¬ 
ages, by using a customized code with an aperture radius 
of 10 pixels (0"20). As a bright star is located at 108 pix¬ 
els (2"2) west of the target and its point-spread function 
(PSF) tail spreads toward the target, we carefully select 
a region to estimate the sky level such that the distance 
from the centroid of the bright star is the same as for 
the target and there is no other noticeable flux contami¬ 
nation. The selected region, a box with the size of 20 x 
20 pixels, is indicated in yellow in Figure [H We calcu¬ 
late a median sky value from this region, and subtract it 
from the fluxes in the target’s aperture. The measured 
target flux is then calibrated to the 2MASS magnitude 
system. For this calibration, we construct a photomet¬ 
ric calibration ladder using J-, H-, and ATg-band archive 
images of VISTA Variables in the Via Lactea (VVV; 
iMinniti et al.ll201^ . because there are only three over¬ 
lapping stars between the IRCS images and the 2MASS 
catalog. The VVV images are three times finer in pixel 
scale and four times deeper in limiting magnitude com¬ 
pared to the 2MASS catalog. We here approximate that 
the IRCS J, H, and K' bands are identical to the 2MASS 
J, H, and bands, respectively. On the VVV images, 
we perform stellar extraction and PSF-fitting photome¬ 
try on a 215 x 215 subarea around the target by using the 
DoPHOT package. Among 40-60 common stars with the 
2MASS catalog, 30 bright-end stars are used for photo¬ 
metric calibration in order to avoid the effect of blending 
on the fainter objects. For the iSMSarM/IRCS images, we 
carefully select 12, 11, and 10 calibration stars for J, 
H, and Kg, respectively, such that they are well isolated 
on the IRCS images and also detected on the VVV im¬ 
ages. Aperture photometry of these stars is done by the 
same manner as for the target object, but this time the 
sky levels are estimated from the stellar-centroid annu¬ 
lus with an inside and outside radii of 70 and 80 pixels 
(ll'4 and 176), respectively. The constructed calibration 
ladders are shown in Figure [6l As a result, we measure 


the target brightness in J, H, and Kg as 


Jtarget = 17.704 ± 0.034 

(11) 

i7target = 17.207 ±0.039 

(12) 

i^s.target = 17.071 ± 0.044. 

(13) 

5.3. The Excess Flux 



The final step to extract the excess flux is to sub¬ 
tract the apparent source fluxes from the measured target 
ones. 

Although the apparent i7-band magnitude of the 
source star is measured from the light curve fitting as 
Hs = 18.569 ± 0.032, no J- or ATg-band light curves 
were obtained. So, we estimate the apparent J- and 
A"s-band source magnitudes Js and Ks,s as follows. 
First, using the co l or-me tallicity-temperature relation of 
iCasagrande et al.l (|2010ll . we obtain the intrinsic source 


colors in (V — J) and {V — Kg) as (V —J)s,o = 1.13±0.05 
and (V — A"s)s,o = 1-48 ± 0.06, respectively. Next, sub¬ 
stituting these values from Vs,o = Hsg — {V — i7)s,o = 
19.73 ± 0.07, we obtain the intrinsic source magnitudes 
of <^S,o = 18.60 ± 0.09 and Kg^s,o = 18.25 ± 0.09, respec¬ 
tively. Finally, the apparent J and Kg source magnitudes 
are derived as Js = 19.02 ±0.09 and Kg^s = 18.40 ±0.10 
by adding Aj and Ak„ to Js,o and ATg^s.O; respectively. 

The best-fit light curve models suggest that the source 
star was still magnified by 1.474 at the time of Subaru 
observation, meaning that the source magnitudes in J, 
H, and Kg at the time of the Subaru observation were 
18.60 ± 0.09, 18.15 ± 0.032, and 17.98 ± 0.10, respec¬ 
tively. All of these magnitudes are significantly fainter 
than the measured target ones, indicating that excess 
flux at the source position is clearly detected. Subtract¬ 
ing the source fluxes from the target ones, we obtain the 
J-, H-, and A'g-band magnitudes of the excess flux: 


= 18.33 ±0.09 

(14) 

= 17.80 ±0.07 

(15) 

= 17.69 ±0.11. 

(16) 


5.4. Origin of the Excess Flux 

In principle, there are three possible scenarios for the 
origin of the detected excess flux: (1) it solely comes 
from the lens star, (2) it is a combination of fluxes from 
the host star and from other astronomical sources (such 
as an unrelated chance-alignment star, a companion to 
the source star, and/or a companion to the lens star), 
and (3) it entirely comes from other astronomical sources 
mentioned previously. 

In the first case, the lens star is almost certainly an 
MS star rather than a giant or stellar remnant, given the 
brightness of the excess flux. In the second case, as we 
will show in Section r6.2l the amount of the contamination 
flux from the extra sources is constrained to < 50% of 
the measured excess flux by the 0 e measurement and the 
upper limit of tte; with this limitation, the lens star is 
also most likely an MS star. 

In the final case, the lens star must not be an MS 
star, but rather a stellar remnant such as a white dwarf 
(WD), a neutron star (NS), or a black hole (BH). This 
could happen in principle, however we ignore this pos¬ 
sibility here because of their relatively low population 
in the Galaxy (M S:WD:NS:BH - 1:0.18:0.021:0.0031; 
iSumi etldl 1201111 and the presumably low planetary 
abundance around them (no planet has yet been detected 
around a WD or BH). In particular, these objects are cre¬ 
ated after disruptive stellar evolutions, including radius 
inflation for all cases and subsequent violent explosions 
for the case of an NS and BH, which could reduce plan¬ 
etary abundance around them. Secondary planet forma¬ 
tion after the evolutions might be possible, however, its 
efficiency is not yet known. 

Therefore, hereafter we simply assume that the lens 
star is an MS star, and consider only the scenarios (1) 
and (2). We note that the stellar remnant scenario can 
be tested by obtaining high-resolution images in the near 
future when the source and lens stars will be separated 
enough; we would not see any object at the expected 
separation if the lens were a remnant. 
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Figure 6. Photometric-calibration ladders for the J (left), H (middle), and Ks (right) bands. For all panels, the black points indicate the 
common stars between IRCS and VVV used to calibrate the instrumental magnitudes of IRCS (bottom axis) to that of VVV (left axis), 
while the red points indicate the common stars between VVV and 2MASS used to calibrate the instrumental VVV magnitudes (top axis) 
to the 2MASS ones (right axis). The location of the target object is indicated as green circle. 


6 . PHYSICAL PARAMETERS OF THE LENS 
SYSTEM 

6.1. No Contamination Assumption 

The absolute physical parameters of the lens system 
can be derived by combining the lens flux and 0 e, both 
of which are independent functions of the host star’s 
mass, TOhost, and the distance to the lens system from 
the Earth, Dp- 

Using 0E, Whost is expressed by 


^host 


^_ 

1 + q KTTrel ’ 


(17) 


where TTrei = AU(1 /L>l — 1/71s)- In Figure [3 we plot the 
nthost"I^L relation for the 0^ value derived in Section lT^ 
fixing Ds at 9.1 kpc. 

The other TOhost~I^L relation can be derived from the 
lens flux. We first assume that the observed excess flux 
comes solely from the lens flux (i.e., there is no contam¬ 
ination from other sources). Because the excess bright¬ 
ness is measured most precisely in the H band among the 
three bands, we use iLexcess to derive the physical param¬ 
eters. The absolute iJ-band magnitude of the host star, 
is calculated using the iJ-band lens flux Hi^, here 
iLp = Hexcess, as a function of Up by the following equa¬ 
tion. 


-M//p = iJp — A//p — 5 log—, (18) 

lOpc 

where Ah.l is the iL-band extinction up to Dp, for 
which we simply assume Ah^l = AnDj^/Dg. For a given 
Up, Mjp^p can be converted to mpost via an isochrone 
model. We use the PARSEC isochro nes version 1.2S 
(jBressan et al.l[2012t iChen et al.l 120141 ). which was im¬ 
proved for low-mass stars over the previous versions. We 
assume the host star’s age of 4 Gyr and a solar metal- 
licity. In Figure [3 we include a plot of this mhost"I^L 
relation. We also plot the 3 ct excluded area from the tte 
upper limit calculated in Section [3^ as the cyan hatched 
region. As a result, we find that the two mhost"I^L re¬ 
lations cross each other at 0.4 < Hp/kpc < 2.6 and 
0.1 < mhost/M 0 < 0.55 within 2 cr, indicating that the 
host star is likely a nearby M dwarf. 

To properly estimate mpost and Up as well as other re¬ 
lated physical parameters including the planetary mass, 


mp, and the projected star-planet separation, rj_, we cal¬ 
culate probability distributions of these parameters by 
means of the Monte Carlo technique. Specifically, we 
solve for mpost, Wp, and r± from a given set of the 
following observed and assumed parameters: X = {p, q, 
s, 0*, Hi^, Ah, Ds, [M/H]}, where [M/H] is the metal- 
licity of the host star. We repeat this calculation about 
10® times by randomly sampling X. For the random dis¬ 
tributions, we use the posterior distr ibut ions obtained 
from the MCMC analysis in Section 13.31 for p, q, and 
s, and we use a Gaussian distribution for 0*, iJp, Ah, 
Ds, and [M/H]. We assume [M/H]=-0.05 ± 0.2, which 
is consistent w ith the metallicity distribution of nearby 
M dwarfs (e.g. iGaidos et aLll2014D . Note that we fix the 
stellar age at 4 Gyr, because the age dependence on the 
isochrone models is negligible for moderately mature M 
dwarfs (>0.5 Gyr). 

The probability distributions of the lens-system pa¬ 
rameters are shown in Figure 13 Each distribution shows 
a bimodal shape, which comes from the fact that the two 
■m-host-I^L relations form a Y-like shape. If 0 e takes an 
upper-side value (above the black solid line), then the 
stellar mass will most likely be very low (^ I. 2 M 0 ), 
which forms the left-hand sharp peak. On the other 
hand, if 0 e takes a lower-side value, then the stellar 
mass can take a value in a wider range of ^0.3-0.5 Mq, 
which forms the right-hand loose peak. Therefore, this 
bimodal shape does not come from two distinct solu¬ 
tions but from one solution. For this reason, for a final 
value of each parameter and its uncertainties, we just 
take the median and 68.3 % confidence intervals of each 
probability distribution. As a result, we obtain ruhost = 
0.34 M 0 , Di^ = 1.3 j kpc, mp = 0.39 Mjup 

(123 1 y 3 MEarth), and r_L = 0.74 I 0 42 AU for the close 
model and 4.3 I 2 5 for the wide model. In addi¬ 
tion, we calculate the probability distribution of the de- 
projected semi-major axis a, assuming a circular orbit, 
uniform distribution of a, a nd random distributio ns of or¬ 
bital inclination and phase (jGould fc Loebl[T992D , result¬ 
ing in the estimated semi-major axis of adrc = O-SO/IqJo 
AU for the close model and Ucirc = 5 . 2 /) 2 'g AU for the 
wide model. All the resultant parameters are listed in 
Table 01 

6.2. Contamination Scenarios 
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Table 4 

Parameters of the Planetary System 


Parameters 

Units 


Values 


Contamination fraction “ 

% 

0 (<10) 

20 (10-30) 

40 (30-50) 

Stellar mass (mhost) 

Mo 

-0.20 

u-iy _o.o7 

O.iO _o.02 

Planetary mass (rrip) 

Mjup 

0.39 

0 22 

-0.08 

n 1 e: -1-0.19 
-0.02 


Me 

123 

70 til 

47 

Distance (Dl) 

kpc 

1 q -|-0.6 
-0.8 

0.85 

0.63 

Projected separation (rj^) 

AU 




Close model 


0 74 

-0.42 

0.45 

r» qq -t-U.dy 

U.OO _o.06 

Wide model 


4.3 tli 

2 ^ +2.Q 
-1.0 

1 q -1-2.3 
-0.3 

Semi-maior axis 

AU 




Close model 


0-90 

n fi7 

^•0/ _o.34 

n 45 

_o.i4 

Wide model 


5-2 ttl 

o q -t-3.2 
-2.0 

2 fi +31 
-0.8 

Contamination probability 

% 




Chance alignment star 


92.9 

5.2 

1.9 

Companion to the source 


89.7 

6.7 

3.5 

Companion to the lens 


99.62 

0.22 

0.16 

Total 


83.2 

11.4 

5.4 


^ The number in parenthesis is for the probability listed in the bottom 
part of this table. 
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Figure 7. Constraints on the distance and mass of the host star. 
The black solid line shows a relation from the best-fit value, 
and the dark gray and light gray regions indicate its la and 2c7 
confidence regions, respectively. The red-bold line is from the 
value with no contamination assumption, and two thin red lines 
are from its 2a upper and lower values. The dashed and dotted 
lines are from the values assuming that 20% and 40% of the 
excess fiux is from contamination sources, respectively. The cyan 
hatched area represents the excluded region calculated from the Sc¬ 
upper limit on tte- 


If the measured excess flux is contaminated by ex- 
tra sources, the actual lens flux would decrease by the 
contaminated flux, pushing down the lens mass. The 
amount of this decrement, however, is limited due to the 
existence of the upper limit of tte. In other words, there 
is an upper limit of the contamination flux; we estimate 
^^50% of the excess flux as this limit above which the two 
mass-distance relations, from Hi^ and from 9^, do not 
cross each other outside the tte excluding region within 


the uncertainties. We look into the effects and probabil¬ 
ities of contaminations in the following sections. 

6.2.1. Effects of Contamination 

First, we estimate how the contamination changes the 
physical parameters of the lens system. We assume two 
contamination levels, / = 0.2 and 0.4, where / is the ratio 
of the contamination flux to the excess flux, and calculate 
the posterior probabilities of the physical parameters for 
each contamination level in the same way as the previous 
section. The results are summarized in Table H) As the 
contamination level increases, the stellar (and planetary) 
mass and distance shrink to ~l/3 and '^1/2, respectively. 
Nevertheless, in many cases the Icr allowed regions of 
each parameter overlap each other; the result obtained 
in Section lOI that the lens system is a nearby M-dwarf 
planetary system, does not change. 

6.2.2. Chance-alignment Star 

Next, we calculate the probabilities of contamination 
for each extra source and for each contamination level. 
The first case of contamination source is an unrelated 
chance-alignment star. On the 18" x 18" FOV of the 
IRCS/iJ-band image, we detect 58 and 53 stars in the 
flux ranges for the contamination levels of / = 0.2 (10%- 
30% of the excess flux) and / = 0.4 (30%-50% of the 
excess flux), respectively. To account for the detection 
incompleteness, we embed a hundred of artihcial stellar 
objects on the IRCS image with a flux corresponding to 
each contamination level, and try to detect them. We 
recover 29% and 72% of the embedded stars for the re¬ 
spective contamination levels, implying that there are 
potentially 197 and 72 stars on the IRCS image, respec¬ 
tively. Adopting a circle with an 8 pixel (0'.'16) radius 
as the contamination cross section, we calculate that the 
chance-alignment probability for the contamination lev¬ 
els of / = 0.2 and 0.4 are 5.2% and 1.9%, respectively. 
The calculated probabilities are summarlized in Table S) 
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Figure 8. Probability distributions of the stellar mass (upper left), distance to the lens system (upper right), planet-star separation for 
the close model (lower left), and that for the wide mode (lower right). In the upper left panel, the planetary mass converted by multiplying 
the stellar mass by g = 1.09 x 10“^ is shown in the top horizontal axis. Probability distributions of the semi-major axis (a) assuming a 
circular orbit, random orbital inclination, and uniform occurrence frequency are indicated in red in the lower panels. In each panel, the 
vertical solid line indicates the median value and gray (red-shared for a) area represents the Icr confidence region. 





6.2.3. Companion to the Source Star 

The second contamination source is a companion to 
the source star. Using Equation (|T^ and the PARSEC 
isochrones version 1.2S, we calculate that 10%-30% and 
30%-50% of the excess flux correspond to the companion 
mass of 0.60-0.83 Mq and 0.83-0.96 Mq, respectively. 
Using the binary (including triplet an d more) fraction 
of the FGK stars of 46% measured by IRaehavan et al.l 
(|2010f) and the companion-mass distribution around the 
FGK stars from the same paper, we estimate that the 
probabilities of a G-dwarf source having a companion 
with the respective mass ranges are 11% and 5.8%. We 
further constrain these probabilities from the limits on 
the projected separation between the source and com¬ 
panion stars. We set an upper limit of 0"16 in the same 
way as the case of the chance-alignment s tar, while set¬ 
ting a lower limit of 1/40 e = 0.36 mas (jBatista et al.l 
1201411 . below which we would see an additional bump in 
the microlensing light curve. This angular range is con¬ 
verted to the semi-major axis range at 9.1 kpc of 4.0 - 
1780 AU, or the orbital-period range of < logP (days) < 
7.4. Adopting a log-normal distribution with the mean 
of logP=5.03 and standard deviation of CTiogp=2.28 as 
the orbital-period dis tribution of FGK-dwarf binaries 
(|Raghavan et al.ll2010ll . we estimate that the fraction of 
binaries within this orbital-period range is 61%. There¬ 
fore, the total probabilities for the contamination levels 
of / = 0.2 and 0.4 are 6.7% and 3.5%, respectively. 

6.2.4. Companion to the Lens Star 


The last case of the contamination source is a compan¬ 
ion to the M-dwarf lens star. In this case, the distance 
to the lens system and the masses of the hypothetical¬ 
binary components change depending on the contamina¬ 
tion level /. We therefore use the following equation to 
calculate the probability that the lens star has a com¬ 
panion: 

-P(/) = Abinary X T! if) ^ (19) 

/ 

where Abinary is the fraction of M dwarfs that have a 
companion, and Fq^ and Fa^ are the fractions of M-dwarf 
companions that have the mass ratio qc and semi-major 
axis Oc, respectively. 

To calculate F g ^, we estimate the companion mass 
using Equation (1181) and an is ochrone model of the 
PHOENIX/AMES-dusty model (jAllard et al.ll200lll for 
a given /. The companion mass varies in the range of 
0.10-0.11 Mq for 0.1 < / < 0.5. Then, the host star’s 
mass for a given / is calculated in the same way as in Sec¬ 
tion [5^23] to derive qg. The qc value varies in the range 
of 0.35-0.94 depending on /. For the Qc distribution, 
we si mply assume a uniform distribution (jJanson et al.l 
I2014f) . We note that because the mass distribution of 
the M-dwarf companion is not yet clear, the choice of this 
simple distribution might be systematics in the probabil¬ 
ity estimation. However, as shown below, the probability 
of contamination from a lens companion is much smaller 
than the other contamination sources, and therefore the 
choice of mass distribution should not affect the total 
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contamination probability. 

To calculate Fa^, we set the limits on the binary sep¬ 
aration. We set the same upper limit of 0"16 as the 
previous cases, which is converted to ac=229 AU for 
/=0.1 and ac=120 AU for /=0.5. For a lower limit, we 
consider a caustic induced by a hypothetical compan¬ 
ion by which additional anomalies could be produced 
in the microlensing light curve. We adopt an upper 
limit on the full width of the hypothetical caustic as 
w = ^qds\ < uq ^ 1.4 X 10“^, where Sc is the projected 
separation of the companion and the host star normal¬ 
ized by 0E- This inequality provides the lower limits of 
the angular separation of 46 mas for / = 0.1 and 75 mas 
for / = 0.5, which correspond to ac=53-65 AU depend¬ 
ing on /. For the distribution of Oc, we adopt a log¬ 
normal distribution with the mean of log U r (AU)= 0.80 
and a standard deviation of criogac=0.48 (j.lanson et al.l 
1201411 . By summing up Equation (IT^ with a step size 
of A/ = 5%, we calculate the probabilities that the lens 
star has a companion as 0.22% and 0.16% for /=10%- 
30% and 30%-50%, respectively. 

6.2.5. Total Probabilities of Contamination 

We calculate that the total contamination probabili¬ 
ties considering the above three scenarios are 11.4% and 
5.4% for / =10%-30% and 30%-50%, respectively. This 
means that the probability that the flux contamination 
fraction is less than 10% and 30% is 83.2% and 94.6%, 
respectively. Therefore, it is most likely that most of the 
excess flux comes from the lens star itself. We note that 
in the case of / < 10%, the Itr allowed ranges of the phys¬ 
ical parameters of the lens system are entirely included 
in those for the non-contamination case. Hereafter we 
simply assume that there is no contamination, and take 
the parameter values calculated in Section rO.ll as our final 
values. Note that the possibilities of the first two con¬ 
tamination sources, the chance-alignment star and com¬ 
panion to the source star, can be tested in the future by 
spatially resolving the source and lens stars. 


6.3. (J — Kg) and (H — K^) colors of the host star 

Because the J- and ATs-band excess fluxes are also mea¬ 
sured, the ( J — Ks) and {H — Kf) colors of the host star 
can be estimated, which can be used as an independent 
check of whether the host star’s mass derived in Section 
IQ is consistent with an M dwarf. We estimate the red- 
de ning for the host s tar by interpolating the online tables 
of iSchultheis et al.l (|2014f) for the Galactic coordinate 
and distance of the lens system, yielding E{J — = 

0.045 ± 0.024 and E{J - = 0.011 ± 0.007, where 

the allowed range of the distance is taken into account. 
Then, assuming that the excess fluxes come solely from 
the host star, we calculate the de-reddening colors as 
(J - iFs)L=0.60 ± 0.14 and [H - iFs)L=0.10 ± 0.13. 
We plot it on the color-color diagram in Figure IHl along 
with the distribution of the metallicity-mea sured nearby 
M dwarfs presen ted inINew ton et a H I12014H and a main- 
sequence track of iBessell fc Brettf i 1988fl . The position of 
the host star is not exactly at the majority of M dwarfs, 
but is more consistent with the region for K dwarfs. How¬ 
ever, considering the large error bars, the measured col- 



Figure 9. [H — Ks) vs (J — Ks) for M dwarfs. The posi¬ 
tion of OGLE-2012-BLG-0563L is indicated by the magenta point 
with error bars. The open circles, triangles, a nd crosses are 
meta llicity-measured nearby M dwarfs presented in IINewton et al.l 
1201411 within the spectral sub-types of M0-M3, M4-M5, and MO- 
MO, respectively. Golors represent the metallicit y ([Fe/H]). The 
black solid line denotes a main-sequence track of iBessell &: Bret3 

iTosa'i . 

ors are still marginally consistent with metal-poor mid- 
to-early M dwarfs. 

We note that the spectral-type estimation from these 
color measurements is more sensitive to systematics than 
the mass measurement from 0 e and Ffp • Although a sys¬ 
tematical change of as large as 0.4 mag in Hp would 
keep the stellar mass within the M-dwarf range, only a 
0.2 mag systematical shift in colors would easily change 
the inferred spectral type. In addition, potential system¬ 
atics in the color measurements are larger than that in 
77 l- Because there is no J- and ATg-band light curves, 
Js and Kg^s are estimated via several calibration pro¬ 
cesses, including the estimation of Aj and Ak^, which 
coul d be a source of systematics in the colors (see Sec¬ 
tion |422]). On the other hand, Hg is directly measured 
from the light curve, which minimizes the uncertainty 
in Hi^. Therefore, it is more likely that the slight dis¬ 
crepancy in the inferred spectral type of the host star 
comes from systematics, or simply statistical errors, in 
the measured colors rather than those in the measured 
mass. We further note that it could also be explained by 
a contamination of a distant early-type dwarf or a disk 
WD, which would cause a bluer shift in the measured col¬ 
ors. However, the probabilities of these scenarios should 
be much lower than those calculated in Section [6. 2. 21 be¬ 
cause the majority of possible chance-alignment stars are 
bulge dwarfs. Instead of identifying the cause of this dis¬ 
crepancy, we will discuss the prospects of improving the 
color estimation in Section [T] 


7. DISCUSSION 

The derived lens parameters indicate that the lens 
system consists of an M dwarf orbited by a Saturn- 
mass planet. Although the planet-star separation is not 
uniquely constrained due to the wide-close degeneracy, 
both solutions locate the planet at cold regions with the 
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Figure 10. Distribution of known exoplanets in the planetary- 
mass vs. stellar-mass plane. The black dots, blue squares, magenta 
triangles, and red circles are planets observed by the radial veloc¬ 
ity, transit, direct imaging, and microlensing methods, respectively. 
The values of microlensing planets are from literature, while those 
of the others are from http://exoplanet.eu. Microlensing planets 
with and without parallax measurements are denoted as open and 
filled circles, respectively, and those for which high-resolution imag¬ 
ing was used to constrain the masses are indicated by large open 
circles. The OGLE-2012-BLG-0563L system is indicated by green. 


estimated semi-major axis close to (for the close model) 
or well beyond (for the wide model) the snow line, which 
we estimate to be AU. Therefore, the planet is lo¬ 
cated in the region where the most microlensing planets 
have been discovered. 

In Figure 1101 we show the distribution of known 
exoplanets in the stellar-mass versus planetary-mass 
plane. OGLE-2012-BLG-0563Lb is the fifth sub- 
Jupiter-mass (0.2 < rUp/Mjup < 1) microlensing 

planet around M dwarfs (0.08 < ruhost ^ 0.55Mq) 
with the mass constrained by either parallax or high- 
re solution imaging, following OGLE-200 6-BLG-109Lb, 
c (|Gandi et al.ll200{^ iBennett et al.l[20int) . OGLE-2011- 
BLG-02 51Lb (iKains et al.ll2013ll . and OGLE-2011-BLG- 
0265Lb ijSkowron et al.ll2014ll . With this new discovery, 
it has become clear that these sub-Jupiter class planets 
form a population around M dwarfs. In contrast, ~l-2 
Jupiter-mass planets are relatively rare around the same 
type of star. This trend is strengthened when consid¬ 
ering the fact that the detection efficiency is higher for 
more massive planets. 

To make this trend more clear, we show the mass distri¬ 
bution of known microlensing planets around M dwarfs in 
Figure rm The upper panel shows a cumulative distribu¬ 
tion function of the observed planets, with a black line for 
all planets and red one for those with a mass constrained 
by parallax and/or high-resolution imaging. The mid¬ 
dle panel shows a histogram for the same sample, where 
the red shaded region indicates Icr (68%) Poisson uncer- 
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Figure 11. Mass distribution of microlensing planets hosted by 
M dwarfs (0.08 < mi,ost/11^© < 0.55). (Top) a cumulative distri¬ 
bution. The thin black line and bold red line are for all planets 
and for the planets with the mass constrained by parallax and/or 
high-resolution imaging, respectively. The location of OGLE-2012- 
BLG-0563Lb is indicated by an arrow. (Middle) the same as the 
top panel, but a histogram. The shaded area indicates the 68% 
poisson confidence region for the histogram of the mass-constrained 
planets. Note that the lower limits shown in the published version 
are erroneous, which are corrected in this figure. (Bottom) the 
same as the middle panel, but the number of planets per bin is 
corrected by a scaled detection efficiency of (g/q’jup)®’®, where q is 
the planet—star mass ratio and Qjup = Mjup/O.35M0. 


tainty for the mass-constrained sample. The lower panel 
shows a “corrected” histogram for which the number of 
planets per one detection is corrected by an approximate 
relative detection efficiency of ((7/gjup)“i where q is the 
planet-star mass ratio and gjup = Mjup/O.SbM©. We 
adopt a = 0.6 for all planets as a mean v alue for the 
centr al-caustic and planetary-caustic events (|Sumi et al.l 
I2010D . In this histogram, the possible paucity of Jupiters 
compared to sub-Jupiters is marginally seen at the 2 cr 
level. We note that simila r trends have also been ob¬ 
served in RV survey s fe.g.. iBonfils et al.ll20I3D and the 
Kepler survey (e.g., iHoward et al.l 1201^ . although the 
orbital regions probed by these surveys are much inner 
than those probed by microlensing. If this trend is true, 
then it could be a consequence of the planetary forma¬ 
tion process in the core-accretion scheme, in which gas- 
accreting planets around M dwarfs can rarely reach a 
Jovia n mass due to the lack of solid and gas materials 
(e.g.. llda fc LinllTOOSi : lAlibert et al.ll2011h . On the other 
hand, several super-Jupiter-mass planets (mp > 2Mjup) 
have also been discovered around M dwarfs by both 
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microlensing and direct imaging. These planets have 
challenged the core-accretion scenarios, and might have 
formed by othe r mechanism s, such as the gravitational 
instability le.g.. lBos'^l2006ll . We note that the steepness 
of planetary-mass function around low-mass stars has al- 
ready been pointed out by seyeral microlensing studies 
(iGould et al.ll201(l iSumi et al.ll201^ : lCassan et al.l[2012l : 

iClanton fc Gaudill2014ll . however, these studies applied a 

simple mass function with a constant slope over all mass 
range (|Gould et al.ll201fll : LSnmi et al.ll201(i iGassan et al.l 
[Ml), or coarsely classified gas planets, treating those 
with a mass ranging 1 < mp/Mj^p < 10 as one group 
([Glanton fc Gaudil 1201411 . We instead demonstrate that 
it has become possible to discuss a hner structure of mass 
function with increasing the number of planet discover¬ 
ies. 

To further clarify the planetary formation process 
around M dwarfs, further increasing the statistics of mi¬ 
crolensing planets in terms of their number and accuracy 
is required. To this end, continuous efforts of not only 
photometric surveys/follow-ups but also high-resolution 
imaging are important. Indeed, OGLE-2012-BLG-0563 
is the first M-dwarf-host planetary event without paral¬ 
lax for which the lens flux is clearly detected, demonstrat¬ 
ing that ground-based near infrared (NIR) AO imaging 
can play a crucial, complementary role to constrain the 
mass of the M-dwarf-host planetary systems. 

Another important aspect of AO imaging in the NIR 
is that it can in principle provide information about the 
metallicity of microlensing M-dwarf host stars from the 
(J — Kg) and {H — Kg) colors. Giant planet-hosting 
sun-like stars tend to be metal rich, which has been 
cited as strong evidence of the core-accretion s cenarios 
le.g.. iSantos et al.l 12004 (Fischer fc Valentill2005[l . Prob¬ 
ing the metallicity of microlensing M-dwarf host stars 
orbited by giant planets can thus be an important tool 
to test the planetary formation models of cold gas giants 
around M dwarfs. As shown in Section 16.31 the mea¬ 
sured [J — Kg) and [H — Kg) colors of OGLE-2012-BLG- 
0563L are marginally consistent with a low metallicity 
mid-to-early M dwarf; if this is true, then this would 
imply that the discovered planet could be a rare sam¬ 
ple being formed by the core-accretion model, or could 
have formed via other formation mechanisms. However, 
the current error bars for these colors are too large to 
be conclusive. These uncertainties mainly come from 
the several calibration steps to derive the source mag¬ 
nitudes, in particular Js and which are not derived 
from light curves but are estimated via the temperature- 
metallicity-color conversions and the extinction laws. As 
described in Section lOl these calibration processes could 
also be a source of systematics. 

These uncertainties can be reduced by future obser¬ 
vations. Because our IRCS/Ji7Ar'-band observation was 
conducted at the time when the source star was still mag¬ 
nified by a factor of 1.47, a second epoch observation (at 
the baseline) with the same instrument/filters will pro¬ 
vide JHK'-h&xiA “light curves”. These light curves can 
directly provide the instrumental (J — K') and {H — K') 
colors of the host star, which can then be converted to 
the colors in the 2MASS system with only color-color cor¬ 
rections. This additional AO observation will provide a 
better estimation of the host star’s colors while avoiding 
the calibration steps for the source star. Ultimately, one 


can obtain a further better estimation of the host star’s 
color by spatially resolving the host star from the source 
star, thus removing the background “noise”. This obser¬ 
vation will be possible in 10 years after the event peak 
with the current ground-based facilities. Note that such 
a spatially resolving imaging will also be able to improve 
the 0E estimation, and hence to refine th e physical pa¬ 
ramet ers of the lens system, as w as done in iBatista et al.l 
(|2015f) and iBennett et al.l (12015f) . This observation will 
be possible in four years, from the event peak, if the Hub¬ 
ble Space Telescope is used. 

8. SUMMARY 

We present the discovery of a microlensing planet 
OGLE-2012-BLG-0563Lb, which was detected 

through intensive photometric observations of a 
high-magnification event. A light curve analysis 

clearly detects the planetary signal of g ~ 10“^ with 
> 1000. On the other hand, we do not detect a 
clear parallax signal in the light curve; we only place an 
upper limit on tte, preventing us from deriving absolute 
physical parameters of the lens system from the light 
curve alone. 

Thanks t o the spectral i nform ation of the source star 
obtained bv iBensbv et al.l (1201311 at a high-magnification 
state, we derive the source’s intrinsic color and magni¬ 
tude with a minimum assumption about the dust extinc¬ 
tion and reddening, and obtain a better estimation of the 
source’s angular radius 0* as well as the Einstein radius 
0E- We also estimate the source star’s distance to be 
Us = 9.11° ® kpc, from the same spectral information. 

To alternatively constrain the physical lens parame¬ 
ters, we conducted a high-resolution JHK'-hacad imag¬ 
ing by using the Subaru/PiO\%S and IRGS instruments 
at the time when the source star was still magnified by 
a factor of 1.47. We successfully detected the excess flux 
from the host star on the source star position, allow¬ 
ing us to derive the absolute physical parameters of the 
lens system by combining it with 0 e estimation and the 
upper limit on tte. We find that the lens system is lo¬ 
cated at 1.3 Ip g kpc from us, and consists of an M dwarf 
(0.34 IqIo ^©) orbited by a Saturn-mass planet (0.39 
-o ’23 Mjup) at the projected separation of 0.74 I®;®® AU 
(close model) or 4.3 AU (wide model). The proba¬ 
bility of contamination in the measured host star’s flux, 
which would reduce the stellar and planetary masses by 
a factor of up to three, is estimated to be 17%. This pos¬ 
sibility can be tested by future high-resolution imaging. 

This is the fifth sub-Jupiter-mass microlensing planet 
around an M dwarf with the mass constrained by par¬ 
allax and/or imaging. The relatively rich harvest of 
sub-Jupiters around M dwarfs is contrasted with a pos¬ 
sible paucity of ^1-2 Jupiter-mass planets around the 
same type of star. This trend could be a consequence 
of the planetary formation process in the core-accretion 
scheme. 

This is also the first M-dwarf-host planetary event 
without parallax for which the lens flux is clearly de¬ 
tected, demonstrating that ground-based AO imaging 
can play a crucial role to constrain the mass of the M- 
dwarf-host planetary systems. In addition, we show that 
NIR AO imaging can, in principle, constrain the metal¬ 
licity of a microlensing M-dwarf host star from (J — Kg) 
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and {H — Kg) colors. Although the current data only 
marginally prefer a low metallicity of OGLE-2012-BLG- 
0563L, further observations will be able to meaningfully 
constrain the metallicity and provide a new insight into 
the formation scenario of this planet. 
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